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ABSTRACT 

We present the rest-frame 200-320 /iin spectrum of the z=3.91 quasar APM 08279+5255, obtained 
with Z-Spec at the Caltech Submillimeter Observatory. In addition to the J = 8— >7 to J =13— >12 CO 
rotational transitions which dominate the CO cooling, we find six transitions of water originating at 
energy levels ranging up to 643 K. Most are first detections at high redshift, and we have confirmed one 
transition with CARMA. The CO cooling is well-described by our XDR model, assuming Li_iookcV ~ 
1 x 10 46 erg s _1 , and that the gas is distributed over a 550-pc sizescale, per the now- favored ^=4 lensing 
model. The total observed cooling in water corresponds to 6.5xl0 9 L Q , comparable to that of CO. 
We compare the water spectrum with that of Mrk 231, finding that the intensity ratios among the 
high-lying lines are similar, but with a total luminosity scaled up by a factor of ~50. Using this scaling, 
we estimate an average water abundance relative to H2 of 1.4xl0 -7 , a good match to the prediction 
of the chemical network in the XDR model. As with Mrk 231, the high- lying water transitions are 
excited radiatively via absorption in the rest-frame far-infrared, and we show that the powerful dust 
continuum in APM 08279+5255 is more than sufficient to pump this massive reservoir of warm water 
vapor. 

Subject headings: galaxies: active — quasars: emission lines — ISM: molecules — instrumentation: 
spectrographs 



1. INTRODUCTION 

1.1. Water as a Molecular Gas Coolant and Radiation 
Field Probe 

Though virtually impossible to observe in the local 
Universe from the ground, H2O is an important con- 
stituent of interstellar gas under the right conditions. 
The species is chemically e nhanced for temperature s 
above a few hundred Kelvin (jDraine fc RobergdllQSlL 
and is sublimated from grain mantles as dust is heated 
to ^100 K. In these conditions, it becomes an impor- 
tant gas-phase repository for oxygen, with abundan ce 
approaching that of CO (e.g. iCernicharo et al.1 (|2006l ) ). 
With water's large dipole moment, the network of ro- 
tational transitions can either be a dominant molec- 
ular coolant (for T > 30 K and rtH, > 10 5 cm~ 3 ; 
Hollenbach fc McKed (fl979l ): lNeufeld fe Melnickl (fl987i ); 
Neufeld fc Kaufman! (| 19931 . and references therein)), or 
a pathway for coupling mid- and far-infrared continuum 
radiation in to the gas if the local con tinuum is sufficiently 
strong (e.g. iScoville fe Kwanl (|1976| ). 

Many of the most important H2O transitions lie in the 
short submillimeter, between 200-300 /jm, which has un- 
til recently been inaccessible in the local universe, and 
has been only sparsely explored at high redshifts via ob- 
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servations in the millimeter band. With the advent of 
broadband spectrometers enabled by high-perfo r mance 
bolometer arra ys (Herschel Sj^E- lGrTffm et al.l (|2010D 
and Z-Spec-e.g. lEarle et al.l (|2006f )). the water spectrum 
can be explored on galaxy scales locally and at high red- 
shift. Water is not an important coolant in the overall 
energy budget of th e gas in the Milky Way (via COBE, 
iFixsen et all (fl999h or even in the nucleus o f the nearby 
starburst galaxy M82 (|Panuzzo et al.l [20Toh . However 
in ex treme ULIRGs s uch as Arp 220 (iRangwala et ail 
1201 lh . and M rk 231 (Ivan der Werf et al.l (120101 here- 
after vdWIO), iGonzalez-Alfonso et al.1 (j2010L hereafter 
GA10)), the most luminous H2O transitions lines can be 
comparable in strength to those of CO. 

Through comparison with water excitation models, 
GA10 show that the water emission lines in Mrk 231 mea- 
sured with SPIRE are not produced by purely collisional 
excitation; the ratio of high-level transitions to the lower- 
level and ground-state transitions is too large, even for 
shock-like conditions. They present a model in which the 
200-300 /im water emission spectrum is pumped by the 
absorption of shorter-wavelength continuum photons, a 
scenario s upported by absorption measu rements with the 
ISO LWS (| Gonzalez- Alfonso et ai1l2008t ). The excitation 
of water is thus coupled closely to the dust continuum in- 
tensity, and GA10 are able to use the ISO + SPIRE lines 
to constrain the size and opacity of the dust continuum 
source in Mrk 231 between 30 and 150 microns, conclud- 
ing that it is ^110-180 pc in size. 

We might expect to find similar processes at work in 
the early Universe, and our wideband 1 mm surveys 
with Z-Spec are revealing H2O in a handful of high- 
redshift objects ( Bradford et al.; 20 09 (hereafter B09), 
ILupu et all 120101 : lOmont et al.l 120111 ). We present here 
a study of APM 08279+5255 at z=3.91, a particularly 
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convenient redshift for characterizing the water emission 
spectrum in the 1 mm band. 

1.2. APM 08279+ 5255 

APM 08279+5255 is an unusual source by any mea- 
sure. It exhibits an apparent total bolometric lu- 
minosity of 7 x 10 15 L Q , dominated by warm AGN- 
heated dust peaking in the rest-frame 15-30 /im regime 
(jRiechers et al.l I2009L and references therein, hereafter 
R09). There is a cold dust component identified in 
the rest-frame far-IR and submillimeter which con- 
tributes ~3% of the total energy. This may be pro- 
duced by star formation, and the [OIII] 88 /xm lumi- 
nosity (1.06x lO 11 /// L©) is consistent with star forma- 
tion contributing up t o ~35% of this far-IR luminosity 
(jFerkinhoff et alJl2010f >. or -1% of the total in the sys- 
tem. APM 08279+5255 has been studied in detail with 
interferometric observations of CO ranging from J = 1— > 
and J = 2— >1 with the VLA (jPapadopoulos et al.l (|200lD . 
more recently R09) to J = 11 — » 10 with the P lateau 
de Bure Interferometer (PdB) ||Weifi et al.l 120071 here- 
after W07). The compact CO morphology (at least at 
J = l— >0) is similar to the image of the optical quasar — 
only just marginally resolved into sub- images at 0.4", and 
R09 model the CO J = 1— ^-emitting gas as a R~550 pc 
disk, magnified by /i=4. We adopt this magnification 
factor, and a luminosity distance of 35,565 Mpc (appro- 
priate for H a = 71, Q M = 0.27, = 0.73). The W07 
CO measurements then translate to an intrinsic molecu- 
lar gas mass of M to tai ~ 1.0 x 10 11 M Q using an adopted 
scale factor of 5.3 x 10 11 M (Kkms -1 pc 2 )" 1 . 

In the R09 lensing model, the magnification is not 
sharply peaked in the central few hundred parsecs: /x(i?) 
ranges from 3-5 out to source radius R ~10 kpc. If 
this is correct, then the CO line ratios and the gas 
conditions inferred from CO (W07) are largely inde- 
pendent of the lensing, and the high-J CO emission in 
APM 08279+5255 is particularly unusual among extra- 
galactic sources. To reproduce the ratios among the low- 
J and high-J CO transitions, W07 require at least two 
components with conditions ranging from cool and dense 
(T~65 K,nn 2 ~lx 10 5 cm~ 3 ) to warmer and less dense 
(T-225 K, n H2 -lxl0 4 cm- 3 ). 

2. OBSERVATIONS 

APM 08279+5255 was observed with Z-Spec at the 
Caltech Submillimeter Observatory (CSO) on a total of 
13 nights among 4 observing campaigns: 2008 March 23- 
28, 2008 April 4-16, 2009 January 6-9 and 2009 Febru- 
ary 23-27. We used the chop-and-node mode with a 20 
second nod period and a 1.6 Hz, 90" chop. The total 
integration time (including chopping) is 25.3 hours, with 
a range of telluric water vapor burdens; T225GHZ ranged 
from 0.04 to 0.20, with a median of 0.083 and a mean of 
0.102. With all data incorporated in a weighted coaddi- 
tion, the sensitivity referred to the total integration time 
ranges from 750 mJy s 1 / 2 in the middle of the band (pro- 
ducing a 2.5 mJy channel RMS) to 1100 mJy sec 1 / 2 at the 
band edges (3.6 mJy channel RMS), with greater degra- 
dation below 200 GHz. The data are calibrated using 
Mars and Uranus, with quasars as secondary calibrators 
using the method described in B09. Figure [T] shows the 
spectrum and the line fits described below. 



2.1. CARMA Follow-Up 

The source was also observed using the 15-element 
Combined Array for Research in M illimeter- wave Astron- 
omy (CARMA, (IBock et al.ll2"006ft h (Figure On 2009 
November 9 and 15, with the array in its C configura- 
tion. The upper sideband was tuned to vjjsb = 236.75 
GHz, corresponding to the H2O 32,i — > 3^2 transition 
at z = 3.9119, the redshift measured with PdB for 
CO J = 9->8 (the highest- J transition measured with 
the PdB) (W07). The lower sideband was centered on 
vlsb = 231.14 GHz. The correlator was configured in 
its 3 x 468.75 MHz setting with partially overlapping 
windows, yielding an effective coverage of Ais sa 1.3 
GHz (Av s» 1600 km s _1 ) with a native resolution of 
fa 31.25 MHz (« 39.6 km s" 1 ). 

The double sideband system temperature was T sys ~ 
280 — 480 K, and the total integration time was 9 hours. 
The quasar 0359 + 509 was the astronomical bandpass 
calibrator (S v ~ 3.5 Jy), and Mars was the absolute flux 
calibrator. The quasars 0927 + 390 (S u ~ 2.3 Jy) and 
0920 + 446 (S„ ~ 1.0 Jy) were primary and secondary 
phase calibrators, respectively, on a calibration cycle of 
10 minutes. The phase RMS on 0927 + 390 after self 
calibration was approximately 25°. Atmospheric decor- 
relation effects were estimated to be 35%, and the fluxes 
have been corrected accordingly. This is not unexpected 
since the observations were obtained in a mid-size con- 
figuration, with typical baselines of ~ 360 m. 

The synthesized beam size is 6 = 0.80" x 0.64" with 
PA = -70°4. Fitting to the LSB and the USB images 
are very consistent and suggest that APM 08279+5255 
is slightly resolved with an intrinsic angular size of 6 « 
0.4", in agre ement with other observations of the source 
angular size (|Krips et al.H2007D . 

3. RESULTS, CONTINUUM AND LINE FLUX 
EXTRACTIONS 

A model spectrum consisting of a a power-law con- 
tinuum with 15 Gaussian spectral lines is fit to the 
160 bolometer flux values, incorporating their measured 
spectral response p rofiles using the method described in 
iNavlor et atl (|2010t ). The redshift is fixed at 3.9119, and 
the line width is Ai>fwhm = 550 km s -1 . This width is 
greater than that reported by W07 for the low-J tran- 
sitions, but is consistent with values measured for the 
J = 9^8 J = 10— >-9 CO transitions. In any case, our 
fits to the Z-Spec spectrum are not very sensitive to the 
adopted line width. The fitted fluxes and resulting un- 
certainties are presented in Table [TJ 

Our CO fluxes are consistent with the W07 IRAM 
measurements to within the uncertainties, though we 
do find systematically larger values, e.g. for J = 9—^8 
and J = 11— >T0 (35-45%). Since the continuum mea- 
surements are consistent to better than this (see Fig- 
ure [!} we rule out global calibration problems at this 
level. It may be that Z-Spec is coupling flux in high- 
velocity wings of the lines which would be missed in the 
IRAM spectra, or there could be spatially extended CO 
emission that the interferometer resolves out. Given the 
coarse resolution, it is also possible that other transi- 
tions near CO are contaminating the Z-Spec fitted flux, 
though aside from the J =10— >9 transition, there are no 
significant candidates in the SPIRE Mrk 231 spectrum 
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Fig. 1. — Z-Spec spectrum of APM 08279+5255. Dashed vertical lines show all transitions of CO and H2O with upper-level energies 
below 700 K which lie in the band. The position of the [Nil] 205 fim transition is also noted. Diamonds show continuum measurements 
from (in increasing frequency): Plateau de Bure (W07), SCUBA ijLewis et alj [l998). CARMA (th is work), CARMA (this work), and the 
SMA (average of two sidebands — horizontal and vertical bars indicate bandwidth and uncertainty) (jKrips et al. 2007). The lower spectrum 
shows the data with a fit to the continuum and CO lines subtracted. 

the two simultaneously, but only their sum is meaningful. 
We estimate the flux of CO J = 10— !>9 as the geometric 
mean of the J = 9— >8 and J = 11— > 10 intensities in 
temperature units, and then subtract this value from the 
fitted flux of the CO+water blend to arrive at an esti- 
mate of the flux for the water line. We find that the water 
and CO fluxes are comparable, though the formal SNR 
for each is modest (~ 2.0) as the errors are propagated 
in this estimate. 

Figure [T]and Table Q] show all eight H2O transitions in 
the band with upper level energies less than 700 K. Of 
these, we detect four in emission with >2.3<7 significance. 
While less than the customary 3a, 2.3a corresponds to 
a likelihood of only (1 - erf (SNR / y/2)) / 2 = 1.0% that 
the transition is not real with positive flux, if the noise is 
Gaussian. In addition to 31,2— >%2,i described above (2a 
is only 2.2% likely to be false), we tentatively identify 
two additional transitions in emission at lower signifi- 
cance (li,i— »0o,o and 2 2i o— >%i,i) and estimate fluxes and 
upper limits. 2o.2- >li,i can only offer an upper limit. 
In Figure 3, we plot 2a upper limits, corresponding to a 
98% confidence level. 

The ground-state fine-structure transition of [Nil] 
( 3 Pi — > 3 Po, A re;jt =205.18 jitm) is identified as tentatively 
in emission at low significance (1.4a). The value is consis - 
tent with the lower limit obtained bv lKrips et al.l ()2007f ) 
with the SMA (3a < 9 Jykms -1 in a single beam beam, 
3a < 16 Jykms -1 if sc aled to the continuum source size 
(|Ferkinhoff et aIll2010h V 
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Fig. 2. — CARMA spectrum of the 32, 1— >3i,2 H 2 transition. 
The velocity scales refers to 2 = 3.9119. The plotted spectra are 
not corrected for the atmospheric decorrelation, though the entries 
in Table 1 are. 



(|van der Werf et all 1201(1 . In any case, our CARMA 
measurement is consistent with Z-Spec in both the line 
flux and continuum to better than the statistical uncer- 
tainties. 

CO J = 10 ^ 9 is clearly blended with the H 2 
3i. 2 — >22,i transition; the two transitions are separated 
by 297 kms" 1 , about 1/3 of a Z-Spec channel. We fit 
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TABLE 1 

1-mm Line Intensity Measurements in APM 08279+5255 
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70.1 
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Z-Spcc fluxes include statistical uncertainties only, IRAM uncertainties arc 


±20%, per W07. 





All luminosities are corrected for lensing, assuming a magnification factor of 4. 
c CO J— 10— 5-9 and H2O 31,2— ^22,1 are blended, see discussion in text. 



the Cloverleaf, and if we assume the same attenuating 
column of Nn, a tt = 3 x 10 23 cm -3 (as is required to 
maintain a molecular phase), then we can use the same 
model as in B09 to consider the CO excitation (Figures 
6, 7, 8 in B09). Assuming the CO emission emerges 
from a 550-pc disk as per R09, then the CO surface 
brightness (ergs" 1 cm" 2 ) in APM 08279+5255 is 0.7- 
0.8 ergs -1 cm -2 , some 3.5-4 times that of the Cloverleaf. 
This emergent flux is readily produced with R<650 pc 
and nH>2xl0 5 cm~ 3 (nn 2 >1 x 10 5 cm -3 if fully molecu- 
lar). While the required density is on the high end of 
that inferred in the W07 CO analysis; this is overall a 
good match to the observed parameters. 

4.2. Water Emission and Radiative Pumping. 

Totaling the available transitions in our band, we find 
a total luminosity of 6.5 xlO 9 L Q in the water rotational 
network. This is a lower bound, but is already close to 
the total CO emission, suggesting that water is energeti- 
cally important to the gas. As GA10 found for Mrk 231, 
the ratio of luminosity in the high-excitation transitions 
relative to the low-lying ones indicates that the high- 
excitation levels are not excited collisionally. The con- 
ditions inferred by W07 using CO (temperature from 
70-200 K and n H2 between 10 4 and 10 5 cm~ 3 ) are well 
below even the extreme conditions considered by GA10 
(T=200 K, iih 2 =1.5x 10 6 cm -3 ) in considering collisional 
excitation. 

Figure [3] shows a comparison of the APM 08279+5255 
H2O spectrum to that of Mrk 231 (adopting 
dL=184 Mpc). To our measurement accuracy, the 
emission in the high-lying lines in APM 08279+5255 
resembles that in Mrk 231, scaled up by ~50x. In 
fact, APM 08279+5255 shows an even higher ratio of 
the E upper > 200 K transitions to the lower-energy 
ones. There is plenty of molecular material in the 
source — the R09 molecular gas mass estimate is a factor 



4. ANALYSIS AND DISCUSSION 

4.1. CO Emission and XDR heating 

We find that the ratios among the high-J CO transi- 
tions are consistent with the excitation model presented 
in W07. If the IRAM fluxes are increased to match the 
Z-Spec values with the factor of 1.4 described above, our 
J = 8^7 and J = 12^11 fluxes and the J = 13^12 
limit all fall on the CO excitation curves in Figures 7 
and 9 of W07, matching to well within the statistical 
uncertainties. We estimate a total lensing-corrected CO 
luminosity of ^7xl0 9 L , a fraction 1.4xl0 -4 of the 
5xl0 13 L output in the far-IR per the decomposition 
of R09. This total CO luminosity fraction is not un- 
usually high; it is comparable to the value in local star- 
burst galaxies and ULIRGs for which the spectra have 
now been measured with Hcrschcl SPIRE. Lco/Lir is 
1.25xl0" 4 in M82 ()Panuzzo et all I2010D . 1.0x10" 



Arp 220 (jRangwala et alj[20TlT) . and 7x 10" 5 in Mrk 231 
(vdWIO). However, as W07 have noted, the CO excita- 
tion of APM 08279+5255 is much higher: the J= 11^10 
to J = 6— t-5 luminosity ratio, for example ranges from 
0.54 ± 0.13 in M82 to - 1 in Mrk 231 to 2.5-3.5 in 
APM 08279+5255, depending on the calibration. Since 
the high-lying CO transitions in Mrk 231 cannot be un- 
derstood in a reasonable PDR framework, the even more 
extreme ratios in APM 08279+5255 cannot be fit by 
PDR models either. 

Given the powerful hard X-ray source, it is likely the 
high-J CO emission i s due to gas cooling in X-ray domi- 
nated regions (XDRs. lMalonev et al.l ()1996[ ) )■ as has been 
modeled for both the Mrk 231 and the Cloverleaf sys- 
tem (B09). For APM 08279+5255, we estimate an in- 
trinsic 1-100 keV luminosity of lx 10 46 ergs -1 , by ex- 
trapolating the apparent luminosity i n the rest-frame 2 - 
10 keV range of 2.3 xlO 46 ergs" 1 per I Just et ail (|2007l ). 
This is similar to the hard X-ray luminosity inferred for 
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of 22 times the GA10 estimate for Mrk 231, including 
all 3 components, but it is some 300 x the modeled 
warm component which dominates the high-lying 
water emission in Mrk 231. Thus even modest water 
abundance throughout, or an abundance similar to that 
obtained in Mrk 231 in 1/6 of the total gas reservoir 
would be sufficient to produce the observed luminosity 
in APM 08279+5255, if it is suitably excited. 

The low significance of the detections and lack of 
higher-frequency coverage for the absorption transitions 
does not warrant a detailed radiative transfer model, but 
it is easy to show that the powerful dust continuum in 
APM 08279+5255 can readily pump the observed water 
transitions with the same mechanism as that modeled by 
GA10 for Mrk 231. Assuming spherical symmetry, the 
observed flux constrains the size of a fiducial region R s 
around a dust emitting source via: 



4< 



d L 



dQ (1) 



where is the luminosity distance, and [i is the magnifi- 
cation factor. J„(f2) is the specific intensity produced by 
the dust distribution; it captures brightness temperature 
and opacity as a function of viewing angle. For a pump- 
ing transition at v em this integral of I v over solid angle 
is the quantity that sets the pumping rate per molecule 
in a given transition. Holding this constant allows an 
estimate of the ratio of the size scale over which similar 
pump conditions exist in APM 08279+5255 as those in 
Mrk 231, through the scaling: 



R s oc d L [{vF u ) obs ^ 



1 1/2 



(2) 



giving 



V w *R s s <xdl [{vF v ) ohs n- v 



3/2 



(3) 

As an example we consider pumping of the 3i_3 level with 
absorption at 58/im, producing the 42,2^4^3 transition 
in emission. Given the observed fluxes, the scaling above 
suggests that R s is 5.1 times larger for APM 08279+5255 
than for Mrk 231. For similar geometries, then, the 
continuum-emitting source in APM 08279+5255 is ca- 
pable of illuminating a volume ~ R? s ~ 130 times larger 
than Mrk 231 with a similar pump rate. This is a rea- 
sonable match to the factor of 50-90 larger luminosity in 
this transition in APM 08279+5255. We expect that for 
APM 08279+5255, the higher color temperature would 
preferentially pump via higher-frequency far-IR transi- 
tions relative to Mrk 231, and this likely explains the 
increased flux in the 52,3— >§i. 4 transition (pumped via 
45 /Ltm) and the reduced flux in 3^2—^22,1 (pumped via 
75 fim). 

The Mrk 231 warm water-emitting component is mod- 
eled by GA10 as a sphere of dust and gas of size 110- 
180 pc. With the above scaling, we know that that simi- 
lar pump conditions must exist in APM 08279+5255 out 
to distances of 560-920 pc. This may seem extreme, 
but APM 08279+5255 is an extremely powerful system; 
consider that the smallest possible physical size for the 
220 K source is in the limit of an optically-thick (at A pea k) 
sphere, for which R m i n ~215 pc. Thus if the molecular 
gas is really situated in a R~550-pc disk per the model 
of R09, then a large fraction of the molecular mass in 
this system is capable of being pumped to the level of 
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Fig. 3. — Comparison of the water spectrum in Mrk 231 with 
that in APM 08279+5255, as measured with Z-Spec. Downward 
arrows show limits at 2<r. The H2O 3i ,2 +2'2, 1 flux is not measured 
in Mrk 231, though as vdWIO note, the J = 10— »9 transition likely 
includes flux from this line. 

the warm component in the Mrk 231 model. Indeed, 
the simple model suggests that the water must exist at 
such large distances because the system is so opaque. 
If we require that the same 220 K, 215-pc dust compo- 
nent account for the emission cit Arest — 250 /xm, then it 
must have an opacity of ^0.8 at this wavelength. Any 
contribution from lower-temperature material requires a 
combination of larger opacity and/or larger physical size. 
In the spherically-symmetric model then, this indicates 
either that the observed fluxes in the water lines must 
be extinction corrected with large factors, or more likely, 
that most of the observed line emission originates at 
R>200 pc. This interpretati on is consistent with that of 
ISrianand fc Petitie an (2000) who find that radiative ex- 
citation from a far-IR source larger than the UV source, 
and at least 200 pc, is required to explain the presence 
of [Sill] and [CII] in the upper far-IR levels indicated in 
optical absorption studies. 

Without more detailed geometric constraints, we can 
only conclude that the total mass of water vapor is at 
least ~50 times that of Mrk 231, some 2.5xl0 4 M , and is 
likely distributed over scales larger than 200 pc, compa- 
rable to the 550-pc size-scale inferred from the CO imag- 
ing. The mass of water vapor corresponds to an average 
abundance relative to H2 of > 1.4 x 10~ 7 , when referred 
to the full molecular gas mass in APM 08279+5255. This 
estimate is a good match to our our XDR model, which 
predicts a water abundance of 1.4-2. Ox 10~ 7 with the 
parameters described above. 
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